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Abstract 

Inclusive D*^ production in two-photon collisions is studied with the L3 detec- 
tor at LEP, using 683 pb~^ of data collected at centre-of-mass energies from 183 
to 209 GeV. Differential cross sections are determined as functions of the trans- 
verse momentum and pseudorapidity of the D*^ mesons in the kinematic region 
1 GeV < Pt < 12 GeV and \r]\ < 1.4. The cross sections a(e+e" e+e"D*=^X) in 
this kinematical region is measured and the a{e'^e~ — >■ e"'"e~ccX) cross section is 
derived. The measurements are compared with next-to-leading order perturbative 
QCD calculations. 
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1 Introduction 



The measurement of open charm production in two-photon colhsions provides a good test of 
perturbative QCD as the large physical scale set by the charm quark mass is expected to 
make the perturbative calculations more reliable. At LEP2 energies, the direct and single 
resolved processes, sketched in Figure [H are predicted P to give comparable contributions to 
the cross section (7(77 —>■ cc). The contributions to open charm production from soft processes 
described by the Vector Dominance Model (VDM) and from double resolved processes are 
expected to be small. The resolved photon cross section is dominated by the photon-gluon 
fusion process 7g — > cc. The cross section for open charm production in two-photon collisions 
therefore depends on the charm quark mass and on the gluonic parton density function of 
the photon. Measurements of open charm production in untagged two-photon collisions were 
performed at PEP 0, PETRA TRISTAN Hj, and LEP [SHE], where charm quarks were 
identified by the presence of D*^ mesons, leptons or Kg mesons. 

This letter describes the results of a measurement of the e^e~ e^e~ccX cross section at 
centre-of-mass energies ranging from 183 to 209 GeV. The data sample corresponds to a total 
integrated luminosity of 683 pb~^ collected with the L3 detector [9^ at LEP. Charm quarks are 
identified by reconstructing D*"*" mesons in three decay channels: 

D*+ ^ (K-7r+)7r+ (1) 

D*+ ^ D°7r+ ^ (K-7r+7r°)7r+ (2) 

D*+ ^ DV ^ (K-7r+7r-7r+)7r+. (3) 

TT^ stands for a "slow pion" , so called as the kinetic energy release in D*"*" decays is only about 
6 MeV. Differential cross sections for D*^ production as a function of its transverse momentum 
Pt and pseudorapidity t] are determined together with the visible cross section. This is then 
extrapolated to derive (T(e+e~ e+e~ccX). 

The PYTHIA ^Uj Monte Carlo program is used to model two-photon processes. It simulates 
77 events according to the current knowledge of hadronic interactions obtained from pp and 
7P studies. Two Monte Carlo samples are used: one with massless quark matrix elements, 
and another with matrix elements with a charm quark mass rric = 1.35 GeV. Both direct and 
resolved diagrams are included in the generation. The resolved processes are generated with the 
SaSld photon parton density function ^1]. A two-photon luminosity function corresponding 
to the Equivalent Photon Approximation J2] is used with a cut off < m^. The L3 detector 
is simulated with the GEANT package ^H] and Monte Carlo events are reconstructed in the 
same way as the data events. Time dependent detector inefficiencies, as monitored during the 
data taking period, are also simulated. 

2 Event Selection and D*"^ Reconstruction 

The event selection follows two steps: first, hadronic final states from two-photon collisions are 
selected, then, D*^ mesons are reconstructed. 

Charge conjugate states are assumed to be included. 
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2.1 Selection of Hadronic Two-Photon Events 



Hadronic events are required to contain at least four particles, with at least three charged 
tracks. A particle is defined as either a track in the central tracker or a calorimetric cluster 
not associated to a track. The background from e"'"e~ — > e"'"e~r+r~ and e"'"e~ r^r" events 
is highly suppressed by this requirement. To reject annihilation events, the visible energy has 
to be less than 80 GeV. The analysis is restricted to events with visible mass VFvis > 2 GeV, as 
calculated from tracks and calorimetric clusters, including those from the small angle luminosity 
monitor. All particles are considered to be pions, except for unmatched electromagnetic clusters 
considered as photons. 

An anti-tagging condition is implemented on the scattered electron or positron, by requiring 
the energy of the most energetic cluster in the forward electromagnetic calorimeter to be less 
than 0.2 ^/s. Only 6% of events in the D*^ signal region are rejected by this criterium. The 
anti-tagging requirement implies a small virtuality of the interacting photons: (Q^) ~ 0.1 GeV^, 
as predicted by the PYTHIA Monte Carlo. 

This preselection yields 4.6 x 10^ events. After applying the above cuts, the contamination 
from annihilation processes is less than 1%. 

2.2 D*^ Reconstruction 

The D*^ selection requires events with less than 15 tracks. Tracks are required to have at 
least 20 hits, a transverse momentum above 0.1 GeV and a distance of closest approach to 
the interaction point in the transverse plane below 1.5 mm. In addition, photon candidates 
are selected in the electromagnetic calorimeter as isolated showers with no close tracks. D*^ 
mesons are selected by first identifying D*^ mesons. Additional criteria reject the combinatorial 
background, and finally, vr^ candidates are added. The D*^ candidates are reconstructed from 
two-track combinations for the channel (^, two-track and two-photon combinations for the 
channel Q, and four-track combinations for the channel Pion or kaon masses are assigned 
to the tracks. We define Pk and as the probabilities for the kaon and pion mass hypotheses, 
calculated from the energy loss, dE/dx, in the tracker. The joint probabilities Pk-Ptt for 
channels (Q) and (j2I) and P^P^P^P^ for channel Q, have to be greater than 0.5% and 
0.1% respectively. In addition, for channel Q the kaon track should satisfy the condition 
Pk/P-it > 2. The vr" candidates for channel Q are constructed from two neutral clusters in the 
barrel electromagnetic calorimeter, with energy greater than 50 MeV and effective mass within 
±20 MeV of the nominal vr*^ mass. The transverse momentum of the 7r° candidate with respect 
to the beam direction must be greater than 400 MeV. 

The effective mass of the D*^ candidate in the K~7r~^, K~7r"'"7r'', or K~7r"'"7r~7r"'" combinations 
is calculated for channels ((H), (j2I) and respectively. This mass should be within ±75 MeV, 
±50 MeV, or ±40 MeV around the nominal D° mass, respectively. These different mass inter- 
vals are determined by the Monte Carlo D*^ mass resolution and reflect the better transverse 
momentum resolution for softer tracks and the high precision vr*^ reconstruction in the BGO 
electromagnetic calorimeter. 

To reduce the combinatorial background in channel (^, the opening angle a between the 
kaon and the D" flight direction has to satisfy the condition cos a > —0.1. 

A kinematical fit using as constraints the D° and tt^ masses is performed for channel Q. The 
fit probability must be greater than 1%. In addition, the sum of the transverse momenta squared 
of these particles must exceed 0.9 GeV . Since the decay D° 7r+7r° proceeds dominantly 

through one of the quasi-two-body intermediate states K*^tt^, K*^7r^ and K~p'^ fT^, the K~7t~^ 
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and K~7r° combinations are required to lie within ±60 MeV of the K*(892) mass and the 7r~7r° 
ones within ±150 MeV of the p mass. This last criterium also holds for channel (jH)), that 
proceeds through the Kvrp state. In addition, for channel Q we require the helicity angle 
9* of an intermediate decay state to satisfy the condition |cos^*| > 0.4, exploiting the P-wave 
properties of vector-particle decay into two scalar particles to reject background. The helicity 
angle 9* is defined as the angle between the direction of a decay product of the vector resonance 
(K*°, K*~ or p+) and the direction of the pseudoscalar particle (vr^, vr^ or K~) from the D° 
decay, calculated in the intermediate resonance rest frame. 

D° decay tracks are required to originate near the interaction point. This is enforced by 
reconstructing a secondary vertex in the transverse plane for channels and Q. The 
impact parameter d with respect to the interaction point in the transverse plane must satisfy 
d < 0.8 mm and d < 0.6 mm respectively, while the D" signed flight distance s should satisfy 
—0.8 < s < 1.0 mm and —0.9 < s < 1.1 mm, respectively. No secondary vertex reconstruction 
is performed for channel but the distance of closest approach of kaon and pion tracks to 
the interaction point in the transverse plane must be less than 1 mm. 

Finally, a slow pion candidate track is added to the combination to form the D*^ candidate. 
The charge must be opposite to that of the kaon candidate, with a distance of closest approach 
to the interaction point in the transverse plane less than 4.5 mm. 

The mass difference AM of the D*^ and D° candidates is shown in Figure El D*^ mesons 
are expected to form a peak around 145.4 MeV, while combinatorial background yields a rising 
distribution, starting at 139.6 MeV, the mass of the charged pion. The reconstructed D*^ 
signal is clearly observed in all the three decay channels under study. The number of observed 
D*^ mesons is estimated by a fit to the AM spectrum. The background shape is modelled by 
Monte Carlo and checked using wrong-charge combinations in data, where the charge of the 
kaon candidate is the same as the charge of the slow pion candidate. The background processes 
e"'"e~ 7/Z qq and e"'"e~ — > W"^W~ — > qqqq represent less than 0.2% of the sample. The 
background shape is modelled by a function with two free parameters a and b: 



where the background shape parameter b is determined from the e+e^ e^e~qq Monte Carlo. 
A fit to the AM spectrum is performed, using a variable width Gaussian to describe the signal, 
and the fitted numbers of D*^ mesons A^^* found in each channel are presented in Table CI 

3 Results 

3.1 Differential cross sections 

The measured differential cross sections for inclusive D*"*^ production as a function of its trans- 
verse momentum Pt and pseudorapidity r], da/dPT and da/d\ri\, are shown in Tables |21 and El 
for the visible region defined by: 



When evaluating da/dP^, which has a steep dependence on P^, a correction factor, determined 
from Monte Carlo, is used to assign the differential cross sections to the centers of the corre- 
sponding bins. The selection efficiency is determined with PYTHIA, using the D*+ branching 
ratios of R eference ITU and assuming a 1 : 1 mixture [T|l7j of direct and single resolved processes. 
The sources of systematic uncertainty are discussed below and summarised in Table |31 



hkgiAM) = aiAM - m^)\ 



(4) 



1 GeV < Pr < 12 GeV and \r]\ < 1.4 . 



(5) 
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• Selection procedure. The uncertainty is estimated by a variation of the selection criteria. 
The largest contributions are the cut on the number of tracks (6.0%) and the D° mass 
window (4.6%) for channel (0); the cut on the transverse momentum of the neutral 
pion (7.8%) and on the sum of transverse momenta (4.3%) for channel the cut on the 
dE/dx probability (11.5%) and the D" mass window (7.5%) for channel Q. 

• Branching ratios. The uncertainties are taken from Reference 

• Background. This uncertainty is obtained by changing the background fit function within 
the uncertainties on its parameters. The range of the fit is also varied. 

• Ratio of the direct to resolved processes. The selection efficiency is re-evaluated changing 
the ratio from 1 : 1 to 3 : 1 or 1 : 3. 

• Monte Carlo statistics. 

• Trigger efficiency. The trigger efficiency and its uncertainty are determined from the data 
using a set of independent triggers. 

Figures El and m show the differential cross sections da/dPT and da/d\T]\. The differential 
cross sections are compared to next-to-leading order (NLO) perturbative QCD calculations, 
based on massive matrix elements In this scheme, the charm quark is not considered 

to be one of the active flavours inside the photon. The Gliick-Reya-Schienbein ^Hl photon 
parton density is used in the calculation. The renormalization scale, fiR, and the factorization 
scale of the photon structure function, jj^p, are taken as fi^ = jj^p/'^ = "^r = a/p^ + m^, 
where pt is the transverse momentum of the charm quark whose mass is rric = 1.5 GeV. The 
calculations are repeated with different renormalization scales separately for the direct and 
single-resolved contributions as well as with different charm quark masses, in order to estimate 
the maximal theory prediction uncertainty. The measurements are in agreement with the NLO 
QCD calculations within the theoretical uncertainties. 

3.2 Visible cross section 

The visible D*^ production cross section a{e~^e~ e~^e~D*^X.)^-^ is calculated in the kinemat- 
ical region ^ as: 

cr(e+e- ^ e+e-D*±X)„, = (6) 

where C is the total integrated luminosity. The trigger efficiency etrig and the selection efficiency 
e^*g are listed in Table Q and the number of observed events in Tables |2l and El Results for 
each Pt and t] bin are shown in Tables El and El while the integrated visible cross sections are 
given in Tabled for each channel. Averaging over the three channels gives: 

a(e+e- ^ e+e"D*^X)^i, = 71.2 ± 5.3 ± 9.8 pb, 

where the first uncertainty is statistical and the second is systematic. 
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3.3 Total cross section 



The total cross section of open charm production is calculated as: 

, + _ + - ^(e^e" e+e-D*±X)^i, - a{e+e- ^ e+e-bbX ^ D*±X)„, 

crfe^e — e^e ccX) = ^ . (7) 

^ ^ ^?o; 2 P(c ^ D*+) ^ ^ 

The probability of a charm quark to hadronize into a D*"*" meson, P(c ^ D*^), is equal to 
0.241 ± 0.008 [17 . The efficiency e^J is a correction factor to extrapolate from the visible 
D*^ kinematical region to the full phase space. It follows from the calculations of Refer- 
ence El as ef^^ = 0.132 

-ao43- The two-photon bb production process also yields D*^ mesons. 
The cross section of this contribution, a{e~^e~ — >• e"'"e~bbX — > D*'^X)ms, is estimated from a 
e"'"e~ e"'"e~bbX Monte Carlo sample and the measured a{e'^e~ e"'"e~bbX) cross section |S] 
as 1.9 ± 0.4 pb. It is subtracted and its uncertainty is included in the systematics. 

The total open charm production cross section at the average centre-of-mass energy of 
197 GeV is estimated to be: 

(T(e+e- ^ e+e-ccX) = (1.12 ± 0.09 ± 0.16 toit ) x Pb- 

The first uncertainty is statistical and second is systematic. The third uncertainty is that on 
the extrapolation from the visible phase space region to the full one. It corresponds to the 
integration [T3j of the dotted and dashed lines in Figures El and El A consistent estimate of 
this uncertainty is also obtained by comparing the extrapolation factors obtained using different 
Monte Carlo generators: the massive and massless PYTHIA samples yield extrapolation factors 
that differ by 24%. 

The total inclusive charm production cross section as a function of the e+e~ centre-of-mass 
energy is shown in Figure El together with our previous measurements jHliHl- The data are 
compared to the NLO QCD predictions of Reference for the direct process and the sum of 
direct and resolved processes. In this calculation, the charm quark mass is taken to be 1.3 GeV 
or 1.7 GeV, and the open charm threshold to be at 3.8 GeV. The resolved process is calculated 
using the Gliick-Reya-Vogt ^Hj photon parton density function. The renormalization and 
factorization scales are taken equal to the charm quark mass. The use of the Drees-Grassie ^H] 
parton density function results in a decrease of the cross section of 9% for rric = 1.3 GeV and of 
3% for nic = 1.7 GeV. Changing the QCD scale from rric to 2mc decreases the predicted cross 
section by 30% for rric = 1.3 GeV and 15% for rric = 1.7 GeV. 

The measured cross section is in good agreement with our previous measurements [OlIHI 
based on lepton tag and with the QCD expectations. The contribution of the single-resolved 
process is needed to describe the data. 
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Channel 


(K-7r+)7r+ 


(K-7r+7r")7r+ 


(K TT+TT 7r+)7r+ 




245 ± 27 
92.9 ± 0.2 
0.61 ± 0.03 


139 ± 18 
92.1 ± 0.3 
0.22 ± 0.02 


99 ± 16 
89.4 ± 0.2 
0.20 ± 0.02 


^(e+e" — >■ e+e^ccX), nb 


63.3 ± 6.7 ± 8.9 
1.00 ± 0.11 ± 0.151^:11 


100.9 ± 13.0 ± 18.2 
1.59 ± 0.21 ± 0.3llg:^^ 


81.1 ± 12.5 ± 17.8 
1.28 ± 0.20 ± 0.29tg:^^ 



Table 1: The number of observed D*^ mesons for the three different channels and the trigger 
{strig) and visible {e^*^ ) efficiencies. Visible and total cross sections are also given. The 
uncertainties on the numbers of reconstructed Y)*^ mesons and on the efficiencies are statistical. 
For the cross section values, the first uncertainty is statistical, the second is systematic. The 
last uncertainty for the total charm cross section corresponds to the extrapolation to the total 
phase space. 



PT(GeV) 




da/dPT (pb/GeV) 


A(T (pb) 


1 - 2 

2- 3 

3- 5 
5-12 


325 ± 31 

130 ± 17 
72 ± 11 
13 ± 5 


54.5 ±5.3 ±6.7 
12.1 ±1.4 ±1.4 
2.9 ±0.4 ±0.4 
0.12 ±0.05 ±0.02 


54.5 ± 5.3 ± 6.7 
12.1 ± 1.4 ± 1.4 
6.7 ± 1.0 ± 1.0 
1.5 ± 0.6 ± 0.3 



Tabic 2: Number of observed D* mesons A'^^^ and measured D* production cross sections 
in different Pt bins for \rj\ < 1.4. Differential da/dPT cross sections refer to the center of the 
bin and Act is the integral over the bin. 





V 






dcr/d ?7 (pb) 


Act (pb) 


0.0 - 0.4 
0.4 - 0.8 
0.8 - 1.4 


194 ± 23 
232 ± 24 
92 ± 17 


49.4 ± 6.0 ± 6.0 

60.3 ± 6.0 ± 7.2 

39.4 ± 7.2 ± 5.6 


19.8 ± 2.4 ± 2.4 
24.1 ± 2.4 ± 2.9 
23.6 ± 4.3 ± 3.4 



Table 3: Number of observed D*^ mesons A^^* and measured D*^ production cross sections 
for 1 GeV < Pt < 12 GeV. Differential d(7/d|77| cross sections refer to the center of the bin and 
Act is the integral over the bin. 
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Source 


(K-7r+)7r+ 


(K-7r+7rO)7r+ 


(K TT^TT 7r^)7rj^ 


Selection 


7.8 


9.9 


16.8 


Branching Ratio [14j 


3.1 


7.2 


4.9 


Fit procedure 


6.6 


6.6 


6.6 


Direct /Resolved ratio 


8.6 


8.6 


8.6 


Monte Carlo statistics 


5.1 


9.4 


8.7 


Trigger efficiency 


0.2 


0.3 


0.2 


Systematics visible 


14.6 


18.9 


22.3 


P(c ^ D*+) 


3.1 


3.1 


3.1 


bb contribution 


0.6 


0.6 


0.6 


Systematics total 


14.9 


19.2 


22.5 



Table 4: Sources of systematic uncertainty in % on the visible D*^ cross section and the total 
charm cross section. The total systematic uncertainty is calculated by adding all contributions 
in quadrature. 



12 




13 



> 





LU 



> 




CO 


-I—' 
LU 



> 




CO 




LU 



150 
100H 
50 




a) 




— backgr.+ signal 
---- backgr. 



L3 Data 



T r 



1 r 



0.14 0.16 0.18 

m(K"7rX) - m(KV) [GeV] 



100 



:r 50- 







b) 



1 r 




backgr.+ signal . ■ o Rata 
backgr. ""^ ^^^^ 



0.14 ^ 0.16 ^ 

nn(K"7cV7Cs) - nn(K"7cV) [GeV] 



0.18 



75 
50 H 
25 




c) 





— backgr.+ signal , ^3 Data 
-- backgr. 



0.14 0.16 0.18 

m(K"jiV7c^7is) - m(K"7C^Ji"ji^) [GeV] 



Figure 2: Distribution of the mass difference D*"*" and D'^ candidates for D*"*" decays into 
a) (K"??"*")??^, b) (K^7r+7r°)7r^ and c) (K~7r"'"7r~7r"'")7r+. The points arc data, the error bars 
represent statistical uncertainties and the hne is the result of the fit to the data points. 
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P-r [GeV] 

Figure 3: The differential cross section da/dPT for inclusive D*^ production. The points are 
the data, the inner error bars represent the statistical uncertainty. The solid curve represents 
the NLO QCD calculations ^S] with the Gliick-Reya-Schienbein jTH] parametrization of the 
parton density of the photon. For the direct (dir) and resolved (res) processes, different values 
of the renormalisation and factorisation scales fiR and hf, as well as of rric, are used. The 
PYTHIA prediction is shown as a histogram. 
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Figure 4: The differential cross section da/dlrjl for inclusive D*^ production. The points are 
the data, the inner error bars represent the statistical uncertainty. The solid curve represents 
the NLO QCD calculations jT^j with the Gliick-Reya-Schienbein jTH] parametrization of the 
parton density of the photon. For the direct (dir) and resolved (res) processes, different values 
of the renormalisation and factorisation scales fiR and hf, as well as of rric, are used. The 
PYTHIA prediction is shown as a histogram. 
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Figure 5: The open charm production cross section in two-photon colhsions as a function of 
the e"'"e~ centre-of-mass energy. L3 measurements based on lepton tag [6,8j are plotted as the 
open circles. The dashed lines correspond to the direct process contribution and the solid lines 
represent the NLO QCD prediction ^Ij for the sum of the direct and resolved processes. The 
effect of different values of rric is also shown. 
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